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Abstract. We construct a multiresolution theory for L 2 (K)©...ffiL 2 (R). For 
a good choice of the dilation and translation operators on these larger spaces, 
it is possible to build singly generated wavelet bases, thus obtaining multires- 
olution super-wavelets. We give a characterization of super-scaling function, 
we analyze the convergence of the cascade algorithms and give examples of 
super-wavelets. Our analysis provides also more insight into the Cohen and 
Lawton condition for the orthogonality of the scaling function in the classical 
case on L 2 (R). 
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1. Introduction 

The applications of wavelet theory to signal processing and image processing are 
now well known. Probably the main reason for the success of the wavelet theory was 
the introduction of the concept of multiresolution analysis (MRA) , which provided 
the right framework to construct orthogonal wavelet bases with good localization 
properties. 

One of the problems in networking is multiplexing, which consists of sending 
multiple signals or streams of information on a carrier at the same time in the form 
of a single, complex signal and then recovering the separate signals at the receiving 
end. 

In HL . Deguang Han and David Larson have shown that the technique of 
multiresolution analysis breaks down, when multiplexing is required, if one just 
amplifies (in the operator theory sense) the steps used in the construction of MRA- 
wavelets. We will state this more precisely in a moment. 

In this paper we will show how multiresolution constructions can be realized for 
multiple signals, provided some slight modifications are done to the usual dilation 
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and translation operators. We believe that our constructions have potential for 
applications to multiplexing problems. 

In |Dut2) and |Dut3| . the second named author introduced a certain affine struc- 
ture on the space L 2 (R) © ... © I? (R) which was shown to admit multiresolution 
wavelet bases. In this paper, we will further analyze this affine structure, give char- 
acterizations of its scaling vectors (Theorem 13.41 Corollary 13.7(1 . Then, the well 
known conditions for the orthogonality of the scaling functions due to A. Cohen 
|Co90| and W. Lawton Law9"Ta] are extended to this larger space in Theorem 13. 91 
We study the convergence of the cascade operator which provides numerical ap- 
proximations for the scaling function. Finally, in Section 5, we construct several 
examples of super-scaling functions and super-wavelets and we prove that all these 
structures admit super-wavelet bases. 

In this introduction we recall several fundamental ideas and notions of the 
wavelet theory. We refer to Dau92 j; |BraJo| for more information on the topic. 

A wavelet is a function ip G L 2 (R) with the property that 

Sjj m T n i) | m, n £ Z} 
is an orthonormal basis for L 2 (R), where U is the dilation operator 

Uf(x) = ±f (|) , (/ei 2 (R),xeR), 

and T is the translation operator 

Tf(x) = f(x-1), (/£L 2 (K),xel). 

A multiresolution analysis is an increasing nest of closed subspaces (V n ) n ^ of 
L 2 (R) which has a dense union, a zero intersection, UV n — V n -i, and there is a 
vector ip in Vo such that 

{T k (f \keZ} 

is an orthonormal basis for Vq. This tp is then called an orthogonal scaling function. 

The wavelets ip are then constructed from the scaling function in such a way 
that 

{T k ip\k S Z} 

is an orthonormal basis for Wq :— V\ © Vq, the relative orthocomplement. 

There are generalizations of the multiresolution analyses some of which we will 
turn to later. One of the directions for these generalizations (see IBCMOp is 
to replace the Hilbert space L 2 (R) by an abstract one H and the dilation and 
translation operators U and T by some abstract unitaries that verify a commutation 
relation such as 

UTIJ- 1 = T 2 . 

The two unitaries will generate what we call a wavelet representation. We adopt 
this representation theoretic view point, and our approach emphasizes the strong 
connections between wavelet theory and the spectral properties of a certain transfer 
operator f Definition II 

We maintain an abstract flavor throughout the paper, but we concentrate on 
some concrete examples of wavelet representations on the Hilbert space L 2 (R) © 
... © L 2 (R) (finite sum), which we will describe in a moment. 

It is known (see |HL| . Proposition 5.16), that no orthogonal scaling function can 
be constructed for L 2 (R)©L 2 (R) with the dilation U®U and the translation TffiT, 
therefore multiplexing can not be obtained in this way. However, we will see that 
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a multircsolution theory can be developed on L 2 (R) ... L 2 (R) with plenty of 
examples of orthogonal scaling functions and wavelets, if some slight modifications 
are done to these operators. 

The dilation and translation operators are constructed as follows: first consider 
a fixed cycle, that is, a periodic orbit for the map z h> z 2 on the unit circle. Let 
C := {zi,Z2,—,z P } be this cycle, zf = Z2,z 2 = z 3 ,...,z 2 _ 1 = z p ,z 2 = z\. The 
Hilbert space is 

H c :=L 2 (M)©...©L 2 (R), 

V v ' 

p times 

the translation operator Tc on He is given by 

Tc(ti, -,C P ) = {ziTti, z p TQ, (6, e i 2 (R)), 
the dilation operator Uc on iJp is defined by 

u c (ti, &) = (u&, u&, cr&, U&), (a, .... e P e L 2 (R)), 

where T and U are the translation and dilation operators on L 2 (R) defined before. 

It was shown in |Dut2l Proposition 2.13] that if a filter m S L°° (T) satisfies 
the condition 

\mo(zi)\ = y/N, (i e {1, ...,p}), 

for some cycle C = {zi, z p } then there is a scaling vector </? S i?c- Such a cycle 
is called an mo-cycle. 

With the operators Uc,Tc, multiresolutions can be constructed and super- 
wavelets (meaning wavelets in spaces larger then L 2 (R)) are obtained. The repre- 
sentations that correspond to different cycles are disjoint so that we will see that 
scaling functions and multiresolutions can be constructed also for the direct sum of 
these representations. Theorem 5.3 shows that, gathering all mo-cycles, one obtains 
an orthogonal scaling vector in the orthogonal sum ®cHc- 

We prove here that the scaling vector constructed out of some m -cycles are 
orthogonal only when all the mo-cycles are taken into consideration, and when this 
is not the case, the super- wavelets form a normalized tight frame (Theorem 13. 911 . 
In particular, one of the consequences of our analysis is that the MRA normalized 
tight frame wavelets, which are known to occur on L 2 (R) when the low-pass filter 
is not judiciously chosen, are in fact projections of good, orthogonal MRA super- 
wavelets. This is due to the fact that the only cycle considered in the construction 
of wavelets on L 2 (R) is the trivial cycle {1}, while the filter mo might have other 
cycles. The direct relation to Cohen's orthogonality condition Co90 is now clear. 

Our paper is structured as follows: 

In Section [5] we define the term of wavelet representation (Definition I2.1fl and 
present the main examples that we will work with (especially Examples 12.51 and 

In Section 13 we define the notions of multiresolution analysis and scaling vector 
(Definition I3.1fl , imitating the one existent for L 2 (R) and give a characterization 
theorem for scaling vectors (Theorem I3.4|) . The theorem is then applied to our 
main examples, and we obtain in this way conditions which characterize scaling 
vectors in L 2 (R) ffi ... L 2 (R) (Corollary l3~7|) . 

Once a scaling vector is found, the construction of wavelets can be done in exactly 
the same fashion as for L 2 (R) (see Proposition 13. 8fl . 
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In Section 3 we will see how one can obtain super-scaling functions aiidsuper- 
wavelets from a trigonometric polynomial filter. Each component of the super- 
scaling function will be compactly supported. Just like in the L 2 (R) case, the 
super-wavelet usually generates a normalized tight frame and, to get an orthogonal 
basis, some extra conditions must be imposed on the initial low-pass filter from 
which the super-wavelet is constructed (such as the Cohen condition or Lawton's 
condition, see Theorem 13.9(1 . 

The scaling vector is approximated by the so-called cascade algorithm. One starts 
with a well chosen function and then the cascade operator is applied successively 
to it. In this way one obtains a sequence which approaches the scaling vector in 
norm. We will see in Section 01 how the initial function must be chosen so that the 
algorithm is convergent. 

Abstract or geometric constructions in wavelet theory must be tested against 
examples and explicit algorithms. We end our paper with Section [S] which con- 
tains several examples showing that plenty of multiresolution super-wavelets can 
be constructed. 

Some notations that we will use in this paper: 
The Fourier transform of / £ L 1 (R) is given by 

/(0 = / f(x)e-**dx. 
Jr 

We denote by T the unit circle {z E C | \z\ = 1} and by /i, the normalized Haar 
measure on T. We often identify functions / on T with 27r-periodic functions on K 
or with functions on the interval [— 7r, 7t]. The identification is given by 

f(z) <-> f{6) where z = e~ ie . 

Many key properties of the scaling vectors are encoded in the spectral properties 
of a certain transfer operator. This is defined as follows: 

Definition 1.1. Let N > 2 be an integer. The transfer operator is associated to a 
function m € L°° (T) and is defined by 

RmJ{z) = ^ E KM| 2 /M, iz € T, / € L 1 (T)). 

W N =Z 

A function h G L 1 (T) is called harmonic with respect to R mo if 

Rm h = h. 

2. Wavelet representations 

Fix an integer N > 2 called the scale. Wavelet representations are an abstract 
version of the situation existent on L 2 (R) as explained in the introduction. The 
Hilbert space L 2 (R) is replaced by an abstract one H and the dilation and trans- 
lation operators are replaced by two unitaries U and T satisfying the following 
commutation relation 

UTU- 1 = T N . 

The translation generates a representation of L°° (T) by the Borel functional cal- 
culus. We give here the proper definition of a wavelet representation and present 
some examples. 
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Definition 2.1. A wavelet representation is a triple tt := (H,U,tt) where H is a 
Hilbert space, U is a unitary on H and 7r is a representation of L°° (T) on H such 
that 

(2.1) UnifjU- 1 = n(f(z N )), (/ € L°° (T)). 

(here, by f(z N ) we mean the map z i— ► f(z N )). 

A wavelet representation is called normal if for any sequence (f n )neM which 
converges pointwise a.e. to a function / e L°° (T) and such that ||/„||oo < Af, 
n e N for some M > 0, the sequence {7r(/„)} converges to tt(/) in the strong 
operator topology. 

Sometimes we call U the dilation and T := tt(z) the translation of the wavelet 
representation (here z indicates the identity function on T, zh>z). 

Example 2.2. The main example of a wavelet representation is the classical one: 
H = L 2 (R), 

tw = ^0 (e^ 9 C*)). 

and 7r is defined by its Fourier transform 

*(/)(0=/£> (/eL°°(T),^L 2 (l)). 

In particular 

T£(x) = e-^^rr), so T(£)(x) = £(x ~ 1), (£ e £ 2 W ,^4 
We denote this normal wavelet representation by SHo- 

Example 2.3. If (Hi, Ui, Hi) are (normal) wavelet representations for i £ {1, n}, 
then (®™ =1 iJj, ©™ =1 C/i, ©™ =1 7Ti) is a (normal) wavelet representation called the direct 
sum of the given wavelet representations. 

Example 2.4. We call cycle a set {zi,...,z p } of distinct points in T, such that 
Zi = z%, Z2 — Z3, Zp = z\. p is called the length of the cycle. {1} is called the 
trivial cycle. 

Let (H,U,tt) be a (normal) wavelet representation. Let C :— {z\, z p } be a 
cycle and a\, a p G T. Define 

H c a ■= H © H e - © H 

S v ' 

p times 

and, for / e (T), e JJ, 

^c,a(Ci. — jCp) : = {aiU&,a2U£,3, ...,a p -iU£p,a p U£i), , 

1"C,a(/)(£l>-,Cp) = ( 7r (/( 2: l 2: ))Cl I 7r (/(2 : 22 : ))6, -,7r(/(%2 ; ))4')- 

Then := (Hc,a,Uc,on^c,a) is a (normal) wavelet representation which we 

call the cyclic amplification of if with cycle C and modulation a. We leave it to 
the reader to check that this is indeed a (normal) wavelet representation (see also 
|Dut2| 1. 

Note that the cyclic amplification with the trivial cycle and ot\ — 1 is the initial 
wavelet representation. 

When ct\ = ... — a p — 1 we will use also the notation ttc '■= 7Tc,q- 
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Example 2.5. If C is a cycle of length p and ai, a p are in T, we denote by 

the cyclic amplification (9to)c,a of the main representation fHo- When all ai are 1, 
we use the shorter notation 

K C = (L 2 (R) a ,U Ci *c). 

Example 2.6. The wavelet representation which is of main importance to us in 
this paper is the direct sum of wavelet representations associated to several cycles. 
That is, if C\, C n are distinct cycles and ai, a n are some finite sets of numbers 
in T, then let C := C\ U ... U C n , a — (oti, a n ) and define 

^c, a ■= OtCx.d © •■■ © 5Kc„,a„, 
which we will call the wavelet representation associated to the cycles Ci, C n and 
the numbers {a;}. 

3. Scaling vectors and compactly supported super-wavelets 

We define now the key concepts of a multiresolution analysis and scaling vectors. 
The definition generalize the existent ones on L 2 (K) (see |Dau92| L 

Definition 3.1. Let (H,U,ir) be a wavelet representation. A sequence (V n ) n£Z of 
closed subspaces of H with the properties 

(i) V n C Vn +X\ 

(ii) U nez y„ = H: 

(iii) n neZ V n = {0}; 

(iv) U(V n ) = V n -H 

(v) There is a <p E V such that {T k ip \ k G Z} is an orthonormal basis for Vq, 
is called a multiresolution analysis (MRA). 

A vector ip 6 H for which there exists a MRA such that (i)-(v) hold, (and <p is 
as in (v)), is called orthogonal scaling vector. 

Before we give a characterization for orthogonal scaling vectors, we need the 
following proposition: 

Definition 3.2. Let 7r = (H, U, tt) be a normal wavelet representation and vi,V2 £ 
H. The Radon-Nikodym derivative h Vl ^ V2 of the linear functional 

/ (tt(/)wi I W2> (/eI°°(T)), 
with respect to the Haar measure on T, is called the correlation function of 
and V2- We also use the notation h v := /it,,^. We can rewrite this as 

(tt(/)«i I «a) = / fh VuV2 dn, (f£L°° (T)). 

The existence of the correlation function is guaranteed by the fact that the 
representation 7r is normal. 

If S is a set of operators on a Hilbert space (such as U plus 7r(/) in the case of 
a wavelet representation), then we denote by S' its commutant (i.e. the set of all 
operators that commute with all operators in S) . 

We proceed towards the theorem that gives a characterization of orthogonal 
scaling vectors. The scaling vector will satisfy a scaling equation ((ii) in the next 
theorem), which relates the scaling vector to a low-pass filter toq £ L°° (T). In 
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order to obtain a MRA, some non-degeneracy conditions must be imposed on mo. 
When mo is degenerate, a residual subspace appears as the intersection of the 
multiresolution subspaces V n (see |EraJo97| 1. 

Definition 3.3. A function mo £ L°° (T) is called degenerate if |mo(^)| = 1 for 
a.e. z€T and there exists a measurable function £ : T — > T and a A £ T such that 

m (z)£{z N ) = X£(z), (z £ T). 

Theorem 3.4. Let tt = (H, U, tt) be a wavelet representation and tp £ H . Then tp 
is an orthogonal scaling vector if and only if the following conditions are satisfied: 

(i) [Orthogonality] The correlation function of tp is h v — 1 a.e.; 

(ii) [Scaling equation] There exists a non- degenerate mo £ L°° (T) such that 
Utp = n(m )ip; 

(hi) [Cyclicity] There is no proper projection p in tt' such that pip = tp. 

Proof. The ideas of the proof are similar to the case of L? (R), so we will only sketch 
them and refer also to |Dau92| . HcWc and BraJo] for more details. 

The orthogonality of the translates is equivalent to the fact that the correlation 
function of tp has Fourier coefficients So, i.e., h v is constant 1. The scaling equation 
is equivalent to the fact that UVq C Vb; the non-degeneracy of mo is equivalent to 
nV^i = {0} (see |BraJo97| and |Jor011 Theorem 5.6]). L)V n is dense if and only if tp 
is cyclic for {U, tt} which is in turn equivalent to condition (iii) 

□ 

Definition 3.5. Let (H, tt, U) be a wavelet representation. A vector tp £ H is called 
a scaling vector with filter mo if it satisfies conditions (ii) and (iii) in Theorem 13. 41 
but here we allow Too to be degenerate. 

Remark 3.6. Note that the condition (iii) is equivalent to 

{U- m ir{f)tp | to > / S L°° (T)} is dense in H, 

since the scaling equation (ii) in 13.41 and the covariance relation 12.1fl are satisfied. 

We will apply the characterization theorem to the instance described in Example 
12.61 When applied to the classical wavelet representation on L 2 (R) we obtain a 
theorem similar to the one in He We , chapter 7. 

Corollary 3.7. Let 9ic,a be the wavelet representation in Examvle \2.b\ Denote 
by e~ l6i ' j the j-th point of the cycle Ci. Then tp c = tp Cl © tp c < 2 © ... © tpc n is 
an orthogonal scaling function for this wavelet representation if and only if the 
following conditions are satisfied: 

(i) 

n Pi 

Y / J2(Peii\p Cl ,k\ 2 ))(^-e i . k ) = 1 for a.e. £ £ R; 
i=l fc=l 

(ii) There exists a function too £ L°° (T) such that for a.e. (£l and for all 
i £ {I, ...,n}: 

a^iVNtfCiANZ) = mo(0i,i + £)fe,i(0. 
a^VNlpCiMN^ = mo(0i,2 + O£c„2(O. 
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(iii) For each i € {1, ...,n}, j G {1, ...,£>$}, fid does not vanish on any subset 
E of R invariant under dilations by N Pi ( i. e. N Pi E — E) of positive 
measure. 

Proof. The formula in (i) is just the correlaton function for h v . Applying the 
Fourier transform to the scaling equation, one obtains (ii). The projections in 
the commutant of this representation are given by sets which are invariant under 
multiplication by N Pi (sec Dut2 , lemma 2.14). The fact that mo is non-degenerate 
is automatic. Indeed, suppose not, then |mo(z)| = 1 for a.e z 6 T so, for all 
ie{l,...,n}, 

VA^|m,i(^ Pl OI = |&7«,l(0l> (£ e «)• 
We will conclude that (pd,i must be 0. If for some a > there is a subset E of 
[— N Pi , —1] U [1, N Pi ] of positive measure such that |£>Ci,i(£)| > a for £ G E, then 

\0 ca (Z)\>-^=, for£eA^£,fceZ 
But then this contradicts the integrability of (pd,i' 

|£c 4 ,i(0| 2 ^>£^iV^A(£) = oc. 

fcGZ 

□ 

As soon as a scaling function is given for a wavelet representation, the construc- 
tion of wavelets follows the procedure described for the L 2 (R)-case in |Dau92) . We 
present below the required ingredients in an abstract version. For a proof look in 
|Dut31 Proposition 5.1]. 

Proposition 3.8. Let n be a normal wavelet representation having an orthogonal 
scaling function <p with non- degenerate filter mo . Denote by (Ki)nez the associated 
MRA. Assume that there are given the "high-pass filters" m\, ..,mN-x £ L°° (T) 
satisfying 
(3.1) 

/ m (z) m (pz) ... m {p N 1 z) \ 
mi(z) mi(pz) ... mi(p N ~ 1 z) 



1 

N 



is unitary for a.e zeT, 

\ m N ^i(z) m N ^i(pz) ... m N -t{p N ~ l z) j 
(p = e^r ), and define ipi E H by 

(3.2) i/> i =:U- 1 ir(m i )'P, (* e {1, .... N - 1}). 
Then 

(3.3) {T fc ^ |* e Z,t € {1,...,JV- 1}} is an orthonormal basis for V\ Vo 
and 

(3.4) {U m T n ipi \m,n £ Z, i e {1,...,JV— 1}} is an orthonormal basis for H. 
Consider a Lipschitz function mo on T that satisfies the following conditions 

(3.5) i? mo l = l, 



(3.6) 



mo has a finite number of zeros. 
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We call a cycle C = {zi = e~ l6>1 , z p — e _j£,p } an mo-cycle if | m-o C-^fe ) I = for 
all k S {1, ...,p}. 

We assume in addition that 

(3.7) There is at least one mo-cycle. 

We have shown in |Dut2| . Proposition 2.13 that for each mo-cycle one can con- 
struct a scaling vector with filter mo in the wavelet representation Die, a where 
Qfe = mo(zk)/yN, (k £ {1, ...,p}). The scaling vector is defined in the Fourier 
space as an infinite product: 

(3.8) Vc,k{x) =[[ ^= -, {x£R),ip c ■= {tpc,i,~.,<Pc,p), 

i=i v N 

(here the subscripts of 9 are considered modulo p, that is 9 a — p , 9_i = 9 p _i, z p+2 = 
Z2, etc.). When m is a trigonometric polynomial each component of the scaling 
vector is compactly supported (see the argument used in Dau92 , lemma 6.2.2). 
As explained in |Dut2j . Proposition 2.13, the function 

p 

h c (6) := h vo {6) = ^Per|£ c , fc | 2 (# - 9 k ), (9 G [-7r,7r]), 
fe=i 

is non-negative, harmonic for R mo and Lipschitz (trigonometric polynomial when 
mo is one). 

Moreover he is constant 1 on the mo-cycle C and it is constant on every other 
mo-cycle. This makes he linearly independent for different m -cycles. 
By remark 5.2.4 in |BraJo| . the dimension of the eigenspace 

{heC(T)\R mo h = h} 

is equal to the number of mo-cycles so that 

{he | C is an mo-cycle} 

is a basis for this eigenspace. This shows that (ii) and (iii) in the next theorem are 
equivalent. Using Theorem 2.16 in |Dut2| we obtain: 

Theorem 3.9. Let mo be a Lipschitz function satisfying \3.5)) . J^S.b}) and and 
let Ci,...,C„ be distinct mo-cycles. Denote by aij the coefficients mo{zij)/\/N 
where for each fixed i the numbers Zij run through the cycle Ci . Define ipc t as m 
13.81) (i G {l,...,n}). Thenipcu '■= fCi ©■•■©<Pc„ is a scaling vector for the wavelet 
representation Vic, a '■= ^Ci,qi © ■■■ © 9^c„,q„ ■ Moreover, the following affirmations 
are equivalent: 

(i) (fc,a is an orthogonal scaling vector; 

(ii) [Cohen's condition] The number of mo- cycles is n; 

(iii) [Lawton's condition] The dimension of the eigenspace 

{h e C(T) | R mo h = h} 

is n. 

If n is smaller then the number of mo-cycles and ipi, ipN-i are defined 
as in Provosition \3.8l then 

(3.9) {WT k ^\te {l,...,N-l},j,keZ}, 
is a normalized tight frame for L 2 (R) . 
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Proof. For the last statement we use the following argument: if n is equal to the 
number of cycles then <p is an orthogonal scaling function so the family in H3.9JI is 
an orthonormal basis. If we project this orthonormal basis onto some of the com- 
ponents corresponding to a choice of a subset of mo-cycles, we get the normalized 
tight frame that corresponds to the case when n is smaller than the number of 
cycles. 

□ 

Remark 3.10. Theorem 18.91 generalizes some well known results of A. Cohen 
and W. Lawton for the classical wavelet representation V\ (see |Co90j . |Law91a| . 
|Dau92| L However much more is true: each normalized tight frame wavelet ob- 
tained in the way described before is in fact a projection of an orthonormal super- 
wavelet, the one which resides in the larger space of the larger wavelet representation 

which takes into consideration all the mo-cycles. 

Of course, it is clear that this projection lies in the commutant of the larger 
representation and so all operators are intertwined. 

4. Convergence of the cascade algorithm 
We saw that the orthogonal scaling vector must satisfy the scaling equation 

Ucp — 7r(mo)iy9. 

Generically, there is no closed formula for the scaling function/ vector. For some 
choices of the filter mo, the scaling function can have a fractal nature (see e.g. 
BraJo ). In applications, numerical values are needed. The cascade algorithm 
provides approximates of the scaling vector ip in the L 2 -norm. It starts with a well 
chosen function V-^°\ an d, by iteration of the refinement (or cascade) operator 

M := U- X -K{m ), 

^O+i) ;= M ^(«), 
it produces a sequence which converges towards the scaling vector tp: 
(4.1) lim \\ip- V? (n) || = 0. 

n — >oo 

The question here is: how should ip(°> be chosen so that the algorithm is conver- 
gent to the scaling vector? We will answer this question in this section. The result 
is related to spectral properties of the transfer operator R mo . For a treatment of 
the L 2 (R) case we refer to |EraJo99] . 

We will consider m , a Lipschitz function on T satisfying (|3.5fl . I|3.tjfl and l|3.7|l . 
Let Ci, ...,C n be all the mo-cycles and define the wavelet representation d\c, a and 
the orthogonal scaling vector ipc,a as is Theorem 13.91 

We will use the notation {zn = e~ t9 ' 1 ,z i 2 = e~ t9z2 , z iPz = e~ l9tp i} for the 
points in the mo-cycle, and the vectors in the Hilbert space of this representation 
are of the form 

with e L 2 (R),i e {l,...,n},j e {1, ...,#}. 
Consider a starting vector in L 2 (l) Co and define inductively 

:= l/c>c,a(mo)V' (n) - 
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The result is 
Theorem 4.1. If 

(4.2) Perl^l 2 ^-^)^^,^,, 
for all £ {1, ...,n},j € {1, ..,?,},/ <E {1, .-,Pi>}, and 

(4.3) $?(2kn) = 5 k , (ie{i,...,n},je{i,..., Pi },keZ), 

then 

(4.4) lim \\ip- V (n) || =0. 

n — >oc 

A simple choice of such a ip^ is 

(4-5) tpij = ~X[o,L), (i € j e {l,...,pi}), 

where L is a positive integer with the property that zf^ = 1 for all i,j. 

Proof. The proof of the theorem will involve the use of some spectral properties of 
the transfer operator R mo regarded as an operator on C(T). 

By Theorem 3.4.4. in BraJo , R mo has a finite number of eigenvalues Ai, X p 
of modulus 1 and R mo has a decomposition 

p 

(4.6) R ma = + S, 

i=l 

where T\ i and S are bounded operators on C(T) such that 

(4-7) Tl = T\ i , T\ i T\. = 0, for i £ j, T Xt S = ST Xi = 0. 

There is a constant M > such that 

(4.8) \\S n \\<M, (neN). 

By Proposition 4.4.4 in Braj^j, the following conditions 14. 9J) and H4.10|) are 
equivalent: 

(4.9) lim RZ (g) = T ^9); 

n—*oo 

(4.10) T Xi (g)=0, for all A, ^1 

By Theorem 2.17 and Corollary 2.18 in |Dut2j . a A with |A| = 1 is an eigenvalue 
for R mo if and only if there is an i € {1, n} such that X Pi — 1. 

For such a A, the eigenspace is described as follows: for each i € {1, ...,n} with 
the property that X Pi — 1, one can define a continuous function with 

Rm Q h^ = Xh^, 

so that for different indices i the functions hf are linearly independent and, in fact, 
they form a basis for the eigenspace 

{heC(T)\R mo h = \h}. 

Define the discrete measures 

(4.11) ^ :=-$> i-1 <W (i e{i, ...,»}, A eT,A« =1), 
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T 



where S z is the Dirac measure at z. Then 

(4-12) T A (/)= ^{f)hf. 

ie{i, ...,«.} with \Pi=i 

These are the tools that we need for the proof of Theorem 14. II 
In the proof we will omit the subscript C, a to simplify the notation. We have 
the following relation between successive approximations: for / E L°° (T) 

/V^+» ^ = (*(f)(<f - ^ (n+1) ) I v - ^ {n+1) ) 

Tr(f(z N )m a )(v - ^ (n) ) 1 n(m )(<p - ^ (n) ) 
/(z JV )|m | 2 /i v _^( n ) 

f{z)R mo h v _^n) dfi. 
Thus, as / is arbitrary: 

This implies by induction that 

= / R mo h tp _^ n -i)d(j t = ...= j Rmoihtp-^m) dp. 

So 

(4-13) ||^-^")|| 2 - / K,. o (V^))^- 

Jt 

Note that, 

h v -Ti,(o) = h v — 2 Re h^^o) + h^o) , 
because, for all real- valued / e L°° (T), 

%{f)(<p - ^ (0) ) | ^ - V (0) ) = WJ> I <p) ~ 2 Re /tt( j> | V (0) \ + U{M {0) \ ^ {0) ) . 



We know that h v — 1 because (/? is an orthogonal scaling vector. 
We see that the correlation function 

n Pi 

(4-14) h^ 0) (6)=Y / J2Pev(!p ij ^f)(d-e ij ), 

i=l j=l 
n Pi 

(4.15) V°> (6) = £ Per ^ | 2 (tf - 

i=i j=i 

We have to compute T\(h v j,(oj) and T\(h^{a)) which amounts to computing the 
values of 7i„ ^,<o) and /i^w on the mo-cycles. For this we will use the inequality 

(4.16) |Per(^gf)(0)| 2 < Per |^-| 2 (0) Per |$°>| a (0). 
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Also, we know (see |Dut2| . Proposition 2.13 and its proof) that the function 
g iS {6) :=Per|^| 2 (6»-%), (i £ {l,...,n},j€ {l,..., Pi }), 

has 

(4.17) 9ij{0kl) = SikSjh 
and also 

(4.18) (p ij {2kn)=5 k (keZ). 
Then, inserting ijHSJ, $M in (EEJ, we get 

(4.19) V°)(^)=Pcr|^ ) | 2 (0) = l. 

Using in (|4.14l) the relations (|4.16|) . 14.17|l and then 14. 3|) and l|4.18|l . we have 



(4-20) h^ m (9 kl ) = Per(^,^ } )(0) = 1. 

With (|4T^|) and tfFZty in (|4~T2|) and JHTTJ, it follows that 

T\(h^( )) = <5a, i • 1 Tx(Re/i ¥3 ^(o)) = <5a,i • 1. 
Therefore (|4.1U|) is satisfied and 

lim #m (V°>) = lim i?" (Reh = 1 

n — too ' ^ n — too ' r ^ 

which implies, with l|4.13|l . that 

lim \\(p-iJj M \\ 2 = 0. 

n — >oo 

It only remains to check that the vector t/>(°J given in equation (|4.5|l satisfies 
(IQl and JO). 



hence 



% 0) (^)=fc, (*€Z). 



Since ^ = 1, we see that e _ii( ^ = 1 so % - 0^-, = ^f- for some fc € Z. 
Therefore, if (i, j) 7^ then fco 7^ mod £ and 



Per IV^m, - 0^) = £ |^)| 2 (^±^) = 0. 

Also 

Per|^| 2 (0) = l 

and 

%?(2kir) = 6 k , (fceZ). 

□ 
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5. Examples 

Example 5.1. Consider the wavelet representation 9\o and suppose mo G L°° (T) 
is a filter that has an orthogonal scaling function <p G L 2 (R) for this representation. 

Define a new filter too G L°° (T) as follows: let p be a positive integer which is 
prime with AT, and define 

m (z) = m (z p ), (zGT). 

Because p and N are mutually prime, the p-th roots of unity, {z G T | z p = 1} split 
into several disjoint cycles (the map z i— > z N is a bijection on {z G T | z p = 1}); for 
example, if N — 2, p = 9 and p^, = e~ l ~, z G {0, ...,8} are the 9-th roots of 1, 
then the cycles are 

{Po}, {pi,P2,Pi,Ps,P7,P5}, {P3, Pel- 
Let Ci, C„ be these cycles 

CiU...UC„ = {zgT|z p = 1}. 
We show that too has an orthogonal scaling vector for the wavelet representation 

namely 

£ := (£ , -,<&>), ^»(ar) = ^ , (x G R, % G {0, - 1}). 

We have to verify the conditions of Corollary 13.71 For this we first write the 
conditions that are satisfied by the scaling function ip: 

(5.1) Per\(p\ 2 (x) = 1, (igR); 

(5.2) \f~N(p{Nx) = m (x)<p(x), (x G R); 



(5.3) There is no ./V-invariant set of positive measure such that (p vanishes on it. 

Having these, we check the conditions for (p. 

The orthogonality condition can be restated in this case as 



^Per|^| 2 (x- — ) = 1 (x£ 
j'=o p 



This holds because 



^Per^l 2 ^ - ^) = \V\ 2 (P* ~ 2J7T + 2knp) = 

^|^| 2 (px + 2fc7r) = Per|£| 2 (px) = 1 [x G R), 

feez 

where we used H5.1(l in the last equality. 

For the scaling equation we have to check that 

N${pNx) = too (p + x^j !p(px), (x G R, j G {0, 1}), 

which is clear from 1)5.2(1 . 
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Suppose now that the cyclicity condition is not satisfied. Then some of the 
components of <p that correspond to one of the cycles do not satisfy the cyclicity 
condition. Let I be the length of this cycle. Then there is an TV'-invariant set of 
positive measure, call it E, such that <p(px) vanishes on E. 

Since (p satisfies the scaling equation, it follows that (p{px) vanishes also on 
NE,N 2 E,...,N l - 1 E. 

Take 

A = -(E UNEU ... U iV^ 1 ^). 
P 

Then A is iV-invariant, of positive measure, and (p(x) vanishes on A. This contra- 
dicts JEU- 

In conclusion, (p is indeed an orthogonal scaling vector with filter rho for the 
wavelet representation 

m Cl ®...®m Cn . 

The next example shows that for the scale N = 2, no matter how we choose the 
cycles C = C\ U C2 U ... U C p , there exists a MRA, orthogonal super- wavelet for the 
representation Dtc- 

Example 5.2. Let d, C 2 , C p be 2-cycles and let C = d U C 2 U ... U C p . We 
will construct an orthogonal scaling vector ip for the wavelet representation 

m Cl ffi ...0 5Hc P - 
The following definitions for i£l will be used in this example : 

(i) x is called a cycle point if there is c in C such that 

x = 9 mod 2-7T, where e~ %e = c; 

(ii) x is called a supplement if x — ir is a cycle point; 

(iii) x is called a main point if it is a cycle point or a supplement; 

(iv) x is called mid-point if x = s ^-, with a, b consecutive main points; 

(v) x is called a cycle midpoint if x — with a, b consecutive cycle points. 

(Here, when we say "consecutive", we refer to the order on the real line ). 
For z e C, z = e-' w \ 9 e [-tt, tt], define 

<Pz{0) — Xr "("o)+»o gq+MM l(fl + ^o)> 

L 2 1 2 J 

where a(^o), #07 b(8o) are consecutive cycle points. It is easy to check that 

Per|£ z | 2 (0-# O ) = 1, (0eR). 

z=e- if> aeC 

Hence ip defined by its Fourier transform 

<fi = ®zeC9z = ©Li ®zed fiz =■ ®i=i<Pd 

is a good candidate for an orthogonal scaling vector corresponding to C build out 
of orthogonal scaling vectors corresponding to each C\ . Next let us define the filter 
m : 

9q£ [— 7r,7r] ,cyclepoint 

where for the cycle point 9q S [— tt, tt], c(6o), 0q, d(6o) are consecutive main points. 
We will now check the scaling equation for the above defined <p and filter mo. It 
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suffices to show that for two consecutive elements zq = e 100 , z\ — e 101 of a cycle Ci 
(i.e. Zq = zi, or equivalently 29q = 9\ mod 2-k) the following holds : 

<p zi (26) = to o (0 + 9 )^(9), a.e. G K 

or, equivalently 

(5.4) <p Zl (29 - 20 o ) = m (9)fi Zo (9 - 9 ), a.e. e R. 

Suppose f> zo (9 — 6q) = Xi a + e o g+go i(^)> where a < 9 < [3 are consecutive cycle 
points. By the definition of too, it follows that there are consecutive main points 
a < 6>o < b such that 

(5.5) m (9)lf Zo (9 - O ) = Xr+ e ° b + e o 1 (6')- 

l 2 ' 2 ) 

Actually a is either a or the first supplement on the left of 9q and b is either (3 or 
the first supplement on the right of 9q. Suppose 5? z , (0 — 6*i) = Xr^iifi h+tn Wi 

l 2 ' 2 I 

with ai < 0i < 6i consecutive cycle points. Then 

l 2 ' 2 J 

and 

(5.6) y, 1 (2fl-2g )= X r .i-»i > 1 -« 1 1 (2g-2e ). 

Note that, under the map x i— > 2x, the consecutive main points a < 9q < b are 
mapped into consecutive cycle points. Since 20o = Q\ mod 27r, and because a 
translation by an integer multiple of 2tt maps consecutive cycle points to consecutive 
cycle points, it follows in particular that there exists fcgZ such that ai — 2a + 2kir, 
9 1 = 20 o + 2kn and b x = 2b + 2kn. We have: 

^i <2e - 2 e < b -^^ 

2 2 

2a + 2kn - (20 o + 2kn) 2b + 2kn - (20 o + 2kn) 

< 20 — 20o < <4> 

2 ~ 2 

a + 9 0<9< b + 9 



2 ~ ~ 2 1 
and, with 1)5. 5|) . I|5.6|l . the relation (|5.4I) is obtained. 

The cyclicity condition is automatically satisfied because all tp z contain a neigh- 
borhood of 0. 

Consequently tp is an orthogonal scaling vector with filter too. 

We can use Example 15.21 to obtain some information about the distribution of 
cycles. The idea is that when Too is the characteristic function of some intervals, 
we can obtain a Cohen condition as in Theorem 13.91 

Proposition 5.3. Suppose too is of the form 

to = s/Nxe, 

where E C [—71", 7r] is a union of intervals such that none of the endpoints of these 
intervals lies on a cycle. Assume moreover that for some distinct cycles Ci, ...,C n , 
the wavelet representation 

has an orthogonal scaling vector with filter mo- 

Then every cycle D disjoint from Ci U ... U C„ must have a point in [—it, tv] \ E. 
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Proof. Since there is an orthogonal scaling vector with filter mo, the condition 

-Rm l = 1 

must be satisfied. 

Suppose there is a cycle D = {e~ Wl , e~ l8p }, different from the given ones 
such that D is contained in E. Because the endpoint of the intervals of E are 
not on cycles, each point of D lies in the interior of E. Then one can construct a 
scaling vector with filter mo for the wavelet representation 9\d, ¥>D,k defined as in 
(|3.8I) . The fact that (pD,k is a well defined L 2 (R) function can be proved using the 
arguments in Dut2 , Proposition 2.13 and Dau92 , lemma 6.2.1. 

The scaling equation can be verified instantly and, since each point of D lies in 
the interior of one of the intervals of mo, it follows that (pD,k is 1 in a neighborhood 
of so the cyclicity condition of Corollary 13. 71 is also clear. 

But then, having a scaling vector with filter mo in the wavelet representation d\r) 
means that 5Hd is the wavelet representation associated to the harmonic function 
h^ D , the correlation function of <po Q.TorOll Theorem 2.4]). 

hg, D k must be bounded for the following reason: clearly <pn,k is either 1 or 0. 
Also, it is impossible to have an x S K with tf>D,k{x) — 1 and tpD^i^lo^) = 1 f° r 
some integer l ^ 0. Indeed, otherwise we would have (from l|3.8|) '): 

mo(^i+B k -i)=VN, m (^j*L + 9 k ^ = Vn, (ieZ,l>l). 

Write lo = N r q with r, q € N, q not divisible by N. Then 

which contradicts i? mo l = 1. 

Thus, if f>D,k(x) = 1 then $>D,k{x + 2kir) = for all integers k ^ 0. This implies 
that Per |^£>,fe| 2 (a;) < 1 so 

h <PD.H ( X ) = £ P(3r \fiD,k\ 2 (x ~ k ) < P . 

Having these, with |Dutll Theorem 2.4], there exists some positive operator 
S that commutes with Uc and nc such that the correlation function hs VClVC 
is h(p D k . Then the correlation function of the vector S x / 2 (pc is h^ D k . Mapping 
C^d" 7I 'd(/) ( ^Z5 to UQ n, Kc{f)S x l' i ipc, w e get a non-trivial operator which intertwines 
the representations associated to D and C. But this is impossible because the 
representations are disjoint |Dut21 Lemma 2.14]. 

In conclusion the assumption was erroneous so the cycle C must intersect the 
complement of E. □ 

We illustrate now some particular cases of the filters, orthogonal scaling vectors, 
and wavelets given in Example 15.21 Here the scale is TV = 2. Since the low-pass 
filter mo is just a characteristic function, the corresponding high-pass filter mi can 
be chosen as the characteristic function of the complement of the set that gives mo, 
and 13.1fl holds. Then the wavelet is defined as in (|3.2I) . 

Example 5.4. Consider the cycle 

^ r -3m. -±zi -6" 

C = {e 7 , e r je r }. 
m = V2x [ _ 7ri _u f]u[ 3 f)W]! 
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2tt 4tt 
7 ! 7 



W=X[-^,f]> ^2 2*], <P3 = X{- 

Then ((pi,tp2, fz) is an orthogonal scaling vector for the wavelet representation 
with filter mo. 

The high-pass filter is 

mi = V2 X[ . 1^,35], 
and the orthogonal wavelet is (^1,^2,^3) with 

^l=X [f ,8-], ^2 = X[_8 f _ f ], "03 = 0. 

Example 5.5. Consider the cycles 

27ri 47ri — 2tti — 4-?ri 2tti — 2-7ri 

Ci = {e 5 , e 5 , e s ?e 5 C 2 = {e 3 , e 3 }. 
m = V^Xb, 

where 

197r n r 15tt IItt, r ll7r 157r n A9tt n 

^ h7r ' "lo" 1 u [ -l3cT' -^T ] u h*' -36-3 u 7r] - 

£?=X[- M , $=Xl-$&Y 
Then ((vi, • ^i), (<Pi , ¥2)) i s an orthogonal scaling vector for the wavelet repre- 
sentation 9tci © ^Hc 2 with filter mo. 
The high-pass filter is 

mi = V^XBi , 

where 

_ 19tt 15tt IItt IItt 15tt 19tt 
and the orthogonal wavelet is ((V>i, V'l); (V*?) "02 )) w ith 

#1=$J=°> V-2 =X[-^,- ft] u[4f,|f] 5 ^4=X[-4 f ,-^]U[ft,^]. 

$=X[ft,ft]. ^ = X[-4 f ,- ft] - 
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